Spatio-Temporal Dependence of the Signaling Response in Immune-Receptor Trafficking Networks Regulated by Cell Density: A Theoretical Model

Cell signaling processes involve receptor trafficking through highly connected networks of interacting components. The binding of surface receptors to their specific ligands is a key factor for the control and triggering of signaling pathways. In most experimental systems, ligand concentration and cell density vary within a wide range of values. Dependence of the signal response on cell density is related with the extracellular volume available per cell. This dependence has previously been studied using non-spatial models which assume that signaling components are well mixed and uniformly distributed in a single compartment. In this paper, a mathematical model that shows the influence exerted by cell density on the spatio-temporal evolution of ligands, cell surface receptors, and intracellular signaling molecules is developed. To this end, partial differential equations were used to model ligand and receptor trafficking dynamics through the different domains of the whole system. This enabled us to analyze several interesting features involved with these systems, namely: a) how the perturbation caused by the signaling response propagates through the system; b) receptor internalization dynamics and how cell density affects the robustness of dose-response curves upon variation of the binding affinity; and c) that enhanced correlations between ligand input and system response are obtained under conditions that result in larger perturbations of the equilibrium . Finally, the results are compared with those obtained by considering that the above components are well mixed in a single compartment

Phosphoproteomics-Based Modeling Defines the Regulatory Mechanism Underlying Aberrant EGFR Signaling

BACKGROUND: Mutation of the epidermal growth factor receptor (EGFR) results in a discordant cell signaling, leading to the development of various diseases. However, the mechanism underlying the alteration of downstream signaling due to such mutation has not yet been completely understood at the system level. Here, we report a phosphoproteomics-based methodology for characterizing the regulatory mechanism underlying aberrant EGFR signaling using computational network modeling. METHODOLOGY/PRINCIPAL FINDINGS: Our phosphoproteomic analysis of the mutation at tyrosine 992 (Y992), one of the multifunctional docking sites of EGFR, revealed network-wide effects of the mutation on EGF signaling in a time-resolved manner. Computational modeling based on the temporal activation profiles enabled us to not only rediscover already-known protein interactions with Y992 and internalization property of mutated EGFR but also further gain model-driven insights into the effect of cellular content and the regulation of EGFR degradation. Our kinetic model also suggested critical reactions facilitating the reconstruction of the diverse effects of the mutation on phosphoproteome dynamics. CONCLUSIONS/SIGNIFICANCE: Our integrative approach provided a mechanistic description of the disorders of mutated EGFR signaling networks, which could facilitate the development of a systematic strategy toward controlling disease-related cell signaling

An empirical Bayesian approach for model-based inference of cellular signaling networks

Background A common challenge in systems biology is to infer mechanistic descriptions of biological process given limited observations of a biological system. Mathematical models are frequently used to represent a belief about the causal relationships among proteins within a signaling network. Bayesian methods provide an attractive framework for inferring the validity of those beliefs in the context of the available data. However, efficient sampling of high-dimensional parameter space and appropriate convergence criteria provide barriers for implementing an empirical Bayesian approach. The objective of this study was to apply an Adaptive Markov chain Monte Carlo technique to a typical study of cellular signaling pathways. Results As an illustrative example, a kinetic model for the early signaling events associated with the epidermal growth factor (EGF) signaling network was calibrated against dynamic measurements observed in primary rat hepatocytes. A convergence criterion, based upon the Gelman-Rubin potential scale reduction factor, was applied to the model predictions. The posterior distributions of the parameters exhibited complicated structure, including significant covariance between specific parameters and a broad range of variance among the parameters. The model predictions, in contrast, were narrowly distributed and were used to identify areas of agreement among a collection of experimental studies. Conclusion In summary, an empirical Bayesian approach was developed for inferring the confidence that one can place in a particular model that describes signal transduction mechanisms and for inferring inconsistencies in experimental measurements

Association Rate Constants of Ras-Effector Interactions Are Evolutionarily Conserved

Evolutionary conservation of protein interaction properties has been shown to be a valuable indication for functional importance. Here we use homology interface modeling of 10 Ras-effector complexes by selecting ortholog proteins from 12 organisms representing the major eukaryotic branches, except plants. We find that with increasing divergence time the sequence similarity decreases with respect to the human protein, but the affinities and association rate constants are conserved as predicted by the protein design algorithm, FoldX. In parallel we have done computer simulations on a minimal network based on Ras-effector interactions, and our results indicate that in the absence of negative feedback, changes in kinetics that result in similar binding constants have strong consequences on network behavior. This, together with the previous results, suggests an important biological role, not only for equilibrium binding constants but also for kinetics in signaling processes involving Ras-effector interactions. Our findings are important to take into consideration in system biology approaches and simulations of biological networks

Signal transduction controls heterogeneous NF-κB dynamics and target gene expression through cytokine-specific refractory states

Cells respond dynamically to pulsatile cytokine stimulation. Here we report that single, or well-spaced pulses of TNFα (>100 min apart) give a high probability of NF-κB activation. However, fewer cells respond to shorter pulse intervals (<100 min) suggesting a heterogeneous refractory state. This refractory state is established in the signal transduction network downstream of TNFR and upstream of IKK, and depends on the level of the NF-κB system negative feedback protein A20. If a second pulse within the refractory phase is IL-1β instead of TNFα, all of the cells respond. This suggests a mechanism by which two cytokines can synergistically activate an inflammatory response. Gene expression analyses show strong correlation between the cellular dynamic response and NF-κB-dependent target gene activation. These data suggest that refractory states in the NF-κB system constitute an inherent design motif of the inflammatory response and we suggest that this may avoid harmful homogenous cellular activation

Systems Biology Modeling Reveals a Possible Mechanism of the Tumor Cell Death upon Oncogene Inactivation in EGFR Addicted Cancers

Despite many evidences supporting the concept of “oncogene addiction” and many hypotheses rationalizing it, there is still a lack of detailed understanding to the precise molecular mechanism underlying oncogene addiction. In this account, we developed a mathematic model of epidermal growth factor receptor (EGFR) associated signaling network, which involves EGFR-driving proliferation/pro-survival signaling pathways Ras/extracellular-signal-regulated kinase (ERK) and phosphoinositol-3 kinase (PI3K)/AKT, and pro-apoptotic signaling pathway apoptosis signal-regulating kinase 1 (ASK1)/p38. In the setting of sustained EGFR activation, the simulation results show a persistent high level of proliferation/pro-survival effectors phospho-ERK and phospho-AKT, and a basal level of pro-apoptotic effector phospho-p38. The potential of p38 activation (apoptotic potential) due to the elevated level of reactive oxygen species (ROS) is largely suppressed by the negative crosstalk between PI3K/AKT and ASK1/p38 pathways. Upon acute EGFR inactivation, the survival signals decay rapidly, followed by a fast increase of the apoptotic signal due to the release of apoptotic potential. Overall, our systems biology modeling together with experimental validations reveals that inhibition of survival signals and concomitant release of apoptotic potential jointly contribute to the tumor cell death following the inhibition of addicted oncogene in EGFR addicted cancers

A multiscale systems perspective on cancer, immunotherapy, and Interleukin-12

Monoclonal antibodies represent some of the most promising molecular targeted immunotherapies. However, understanding mechanisms by which tumors evade elimination by the immune system of the host presents a significant challenge for developing effective cancer immunotherapies. The interaction of cancer cells with the host is a complex process that is distributed across a variety of time and length scales. The time scales range from the dynamics of protein refolding (i.e., microseconds) to the dynamics of disease progression (i.e., years). The length scales span the farthest reaches of the human body (i.e., meters) down to the range of molecular interactions (i.e., nanometers). Limited ranges of time and length scales are used experimentally to observe and quantify changes in physiology due to cancer. Translating knowledge obtained from the limited scales observed experimentally to predict patient response is an essential prerequisite for the rational design of cancer immunotherapies that improve clinical outcomes. In studying multiscale systems, engineers use systems analysis and design to identify important components in a complex system and to test conceptual understanding of the integrated system behavior using simulation. The objective of this review is to summarize interactions between the tumor and cell-mediated immunity from a multiscale perspective. Interleukin-12 and its role in coordinating antibody-dependent cell-mediated cytotoxicity is used illustrate the different time and length scale that underpin cancer immunoediting. An underlying theme in this review is the potential role that simulation can play in translating knowledge across scales